clc;clear;

V_Line = 220e3;

X = 8;
C = 0.025e-6;

R = 600;

f = 50;
w = 2*pi*f;

L = X/w;

% Part 1
w0 = 1/sqrt(L*C);
f0 = w0/2/pi;

% Part 2
alpha = 1/(2*R*C);
fd = (1/2/pi)*sqrt(w0^2-alpha^2);

% Part 3
R_Critical = 0.5*sqrt(L/C);

% Part 4
syms R
Equ = 0.25*f0 == (1/2/pi)*sqrt((1/L/C)-(0.5/R/C)^2);
R = eval(solve(Equ,R));



